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Abstract: We perform a detailed analysis of the detection prospects of neutralino dark 
matter in the mSUGRA framework. We focus on models with a thermal relic density, 
estimated with high accuracy using the DarkSUSY package, in the range favored by current 
precision cosmological measurements. Direct and indirect detection rates are computed 
implementing two models for the dark matter halo, tracing opposite regimes for the phase 
of baryon infall, with fully consistent density profiles and velocity distribution functions. 
This has allowed, for the first time, a fully consistent comparison between direct and indi- 
rect detection prospects. We discuss all relevant regimes in the mSUGRA parameter space, 
underlining relevant effects, and providing the basis for extending the discussion to alter- 
native frameworks. In general, we find that direct detection and searches for antideuterons 
in the cosmic rays seems to be the most promising ways to search for neutralinos in these 
scenarios. 
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1. Introduction 

The identification of dark matter in the Universe is one of the most compelling targets 
in Science today. In the "concordance" cosmological model |ffl], emerging from precision 
cosmological measurements and from tests of the theory of structure formation, some un- 
known form of non-baryonic cold dark matter (CDM) accounts for about 30% of the mean 
energy density of the Universe today. Among the solutions of the dark matter puzzle, the 
most natural scheme is the one in which CDM, analogously to the baryonic and radiation 
components, appears as a thermal relic from the early Universe; in particular, weakly inter- 
acting massive particles (WIMPs) are natural candidates, as their thermal relic abundance 
is automatically of the right order of magnitude. Moreover, it is a scheme with strong 
motivations from the particle physics point of view: the most widely studied WIMP dark 
matter candidate is the lightest neutralino in supersymmetric extensions of the Standard 
Model of particle physics. 

Several techniques have been studied to search for dark matter WIMPs (for thorough 
reviews and comprehensive lists of references, see, e.g. (2|, [||]). One of the issues that is 
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often raised, eventually to understand in what direction experimental efforts should be 
focussed, regards the comparison between capabilities of different techniques. Actually, 
as we already mentioned, the idea of WIMPs is associated to a scheme rather than to a 
model, and no simple recipe can be given. It is only when the setup of a specific model 
is fully defined that a sensible comparison can be performed; on general grounds, the 
only statement one can safely formulate is that, in most cases, different techniques probe 
different properties of WIMP dark matter (we refer here to properties related both to 
particle physics and astrophysics), hence they are complementary. On the other hand, it is 
true that focussing on a specific model is sometimes a very useful exercise, which allows for 
better understanding of at least some aspects of an otherwise too complicated problem. In 
the context of supersymmetric extensions to the Standard Model, this is often done, e.g., 
to present limits from current sets of data, or to address the reach of future accelerator 
experiments; the same is useful applied to dark matter searches. 

In this paper we will focus on the supersymmetric neutralino as a dark matter candi- 
date in the so-called mSUGRA framework of supersymmetry breaking, one of the simplest 
and most popular models. Direct and indirect detection of dark matter in this model has 
been the subject of several studies (a non-exhaustive list of recent papers on this topic 
includes, e.g., @, ||, ||, @, ^, ||, |ll]]). Our reiteration starts from the point of view of 
restricting to configurations corresponding to relic abundances in the range currently fa- 
vored by cosmological measurements; to fulfill this requirement, we calculate relic densities 
with the DarkSUSY computer code [p^| , which includes all possible sfermion, neutralino 
and chargino coannihilations and is currently the most accurate code for relic density cal- 
culations. We will then consider two limiting cases of profiles for the dark matter halo of 
the Milky Way, consistent with available constraints, and derive self-consistent density and 
velocity distributions for the neutralinos in two such sample halo models. This approach 
will allow, for the first time, a fully consistent comparison of the prospects to detect relic 
neutralinos with direct detection or various indirect detection methods, like neutrinos from 
the Earth/Sun or cosmic rays from annihilations in the galactic halo. Our aim is also to 
try to make this comparison more transparent than in previous studies, and to underline 
step by step what the properties are that enter critically to make the balance bend on one 
side or another. This analysis will then try to clarify what the characteristics are of the 
model we are considering, and at the same time to help in foreseeing what may change in 
other setups. 

The outline of the paper is as follows. We will start with an introduction to the 
supersymmetric framework, continue with a discussion of the self-consistent halo profiles 
we use, go through the various detection rates for models within the measured relic density 
range and finally end with a discussion and conclusions. 



2. The particle physics model 
2.1 The supersymmetric setup 



The present analysis is performed in the mSUGRA setup [13], namely the N=l supersym- 
metric extension of the Standard Model (SM) defined in a supergravity inspired framework 
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and with the smallest possible number of free parameters: on top of a structure with min- 
imum field content, universality is assumed at the grand unification (GUT) scale, both in 
the gaugino and the scalar sector of the theory. The mSUGRA action is then fully defined, 
by only four parameters and one sign: the GUT scale values of the soft supersymmetry 
breaking fermionic mass parameter m^, the scalar mass parameter mo and the trilinear 
scalar coupling Aq, the ratio of the vacuum expectation values of the two neutral compo- 
nents of the SU(2) Higgs doublets tan (3, and the sign of the Higgs superfield parameter /i 
(the absolute value of \x is fixed by electroweak symmetry breaking; regarding the conven- 



tion on its sign, following, e.g., |l4j| , it is assumed here that \x appears with a minus sign 
in the superpotential). 

The appropriate set of renormalization group equations (RGEs) allows to relate univo- 
cally the GUT scale structure to the low energy (weak scale) spectrum of the theory. Here, 
soft breaking parameters, gauge and Yukawa couplings are evolved down to the weak scale 
with the ISASUGRA RGE code as given in version 7.67 of the ISAJET software pack- 
age (introducing some minor changes, such as, the conversion of the ISASUGRA code 
to double precision to improve on its stability; for more details on this and other technical 
points on the code implementation and on the interface with the DarkSUSY package, see 



the discussion in [16|). The mSUGRA setup is probably the most popular framework for 
studying supersymmetric (SUSY) extensions of the Standard Model, and its rather con- 
strained low energy structure has been extensively discussed. We review here very briefly 
those features which will be relevant in our discussion of the detection prospects for SUSY 
dark matter in this framework. 

Our working hypothesis is that the lightest supersymmetric particle (LSP) (and our 
dark matter candidate) is the lightest neutralino, defined as the lightest mass eigenstate 
from the superposition of the two neutral gaugino and the two neutral Higgsino fields: 

X? = N XX B + A 12 TU 3 + N 13 H% + A 14 F 2 ° . (2.1) 

The coefficients N\j, obtained by diagonalizing the neutralino mass matrix, are mainly a 
function of the bino and the wino mass parameters M\ and M2, and of the parameter (i. 
From the assumption of gaugino mass unification at the GUT scale, it follows that the 
hierarchy between M\ and M2 at the weak scale is fixed to about M\ ~ O.5M2, and hence 
that the wino component of the LSP is always very small. There are then essentially just 
two regimes in the composition of the lightest neutralino: in most of the allowed regions 
in a generic scan in the mo-m^j parameter space, the neutralino LSP is a very pure bino; 
in a rather thin slice at tuq 3> mi/2, sometimes dubbed the "focus point" region [17], on 



the border with the region where there is no radiative electro-weak symmetry breaking, 
the parameter [i is driven to small values and forces a mixing between the gaugino and 
Higgsino sectors, with a considerable Higgsino fraction in the lightest neutralino. 

Regarding SUSY scalars, except at very large values of mo, i.e. for very heavy sfermions 
which do not play much of a role, the RGEs drive the slepton sector to be lighter than the 
squark sector: in particular in the cosmologically interesting regime mo < m^, the lightest 
stau is always the lightest sfermion, possibly even lighter than the lightest neutralino if 
mo <C m 1( / 2 ; selectrons and smuons are slightly heavier, while the lightest stop, which is 
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also the lightest squark, is significantly heavier and scalar partners of the light quarks, 
whose exchange can contribute to the scattering of neutralinos on protons and neutrons, 
are even more massive. A notable exception to the pattern just described is given by the 
case when the trilinear coupling, A, in the stop mass matrix is tuned to such a value that 
one of the two stop mass eigenstates becomes light: in this case, a stop can become the 
next-to-lightest SUSY particle or even the LSP, while the spectrum of other SUSY scalars 
is not significantly changed. 

Finally, regarding the Higgs sector, the mass of the light CP-even Higgs boson H®, 
especially for low tan/3 and low mo, is mainly set by m-i/2 and changes only slightly 
with mo, while the masses of the heavy CP-even Higgs boson Hq and of the CP-odd Hq 
increases monotonically with both mo and m^. In most regions of the parameter space 
Hq decouples and behaves like the Standard Model Higgs; at large tan (3, the mass of the 
heavy Higgs bosons can be driven to rather small values, and in a diagonal stripe of the 
parameter space m^-mo) the mass of the LSP happens to match one half of the mass of 
the heavy Higgs boson(s), giving rise to resonant annihilations of neutralinos. 

We will focus on a few selected scans in the mSUGRA parameter space, encompassing 
however all these features, whose role on dark matter detection we will enlighten. 

2.2 Constraints on the model 

The first pattern of model discrimination we apply is based on the calculation of the 
relic abundance of the neutralino LSP for any given set of the parameters, and by the 
requirement that the relic density matches the best fit value from the latest cosmological 
measurements of the non-baryonic dark matter component in the Universe. 

The relic density calculation we implement for this analysis is included in a new ex- 
tended version of the DarkSUSY package [l^], which has been recently released. It is suitable 
for a proper treatment of any coannihilation effect, applying, at the same time, the state 
of the art technique to trace the freeze-out of a species in the early Universe | |18| |, with a 
careful numerical treatment of resonance and threshold effects, and full numerical solution 
of the density evolution equation (avoiding approximations such as, e.g., the expansions in 
powers of the relative velocity that is often applied). This method is described in detail in 



Ref. [16], where it has been applied to the SUSY setup considered here; assuming masses, 
widths and couplings of particles in the model are given with an adequate precision, this 
neutralino relic abundance calculation has an estimated precision of 1% or better. 

The combined analysis of the latest data on the cosmic microwave background ani- 
sotropics and large scale galaxy surveys gives a fairly accurate estimate for the cold dark 

matter contribution to the energy density of the Universe. We take, as reference value, 

2 _ n mQ+0.020 

-0.022 



the best fit value derived in Ref. [19], under standard assumptions: QcDAfh 2 = 0.103 



(Table 3, first column, in |19|), where Sic dm is the ratio between mean CDM density and 
critical density and h is the Hubble constant in units of 100 km s _1 Mpc -1 . 

So far, experimental searches for SUSY particles have given null results. To exclude 
models that violate accelerator constraints, we adopt the compilation of limits on SUSY 
masses as given by the Particle Data Group 2002 (PDG) ^(|. On top of these, we consider 
two more restrictive conditions: We assume as limit on the lightest chargino mass the final 
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kinematical limit of LEP, i.e. 103.5 GeV [21]; also, as we mentioned that H® is in most 
case SM-like, we will use as a limit on the i?^ 1 mass, the current limit on the SM Higgs 
boson mass, 114.1 GeV p2fl . 

Another strong bound is given by the SUSY contributions to the rare decay b — ► 57. 
Our estimate of this process includes the complete next-to-leading order (NLO) correction 
for the SM contribution and the dominant NLO corrections for the SUSY term. The NLO 
QCD SM calculation is performed following the analysis in Ref. [23], modified according 
to [ ]24| , and gives a branching ratio BR[1? — > X s j] = 3.72 x 10~ 4 for a photon energy 
greater than mfe/20. In the SUSY contribution, we include the NLO contributions in the 
two Higgs doublet model, following |pfif| , and the corrections due to SUSY particles. The 
latter are calculated under the assumption of minimal flavour violation, with the dominant 
LO contributions from Ref. [26], and with the NLO QCD term with expressions of [27] 
modified in the large tan/? regime according to [^qj . In the mSUGRA framework (see, e.g., 
f28[| ), the largest discrepancy between the LO and the NLO SUSY corrections are found for 
sign// > 0, large tan/? and low values of m^: in this case the SUSY contribution to the 
decay rate is negative, and the discrimination of models based on the NLO analysis is less 
restrictive than the one in the LO analysis. We will assume as allowed range of branching 
ratios 2.0 x 10~ 4 < BR[£? — >• X s 7] < 4.6 x 10~ 4 , which is obtained adding a theoretical 
uncertainty of ±0.5 x 10 -4 to the experimental value quoted by the Particle Data Group 
2002 |0|. 

Finally, we consider the information on the model following from the SUSY contri- 
bution to the anomalous magnetic moment of the muon. Our estimate of the (g — 2)^ 
SUSY term is based on the analysis in Ref. 2£]. This can be compared with the range of 
values favored by the latest BNL measurements, relative to the calculation in the standard 
model H], -5 x 10~ 10 < a* p - af* < 47 x lO" 10 [2a]. The Standard Model prediction 
depends on experimentally measured quantities that can be obtained from r decay or e + e~ 
annihilation data. As these two methods give different results, we have, to be conservative, 
taken the lower limit derived from r data and the upper limit from e + e~ data. This gives 
the rather large range of allowed {g — 2) At -values stated above. With this broad range of 
allowed values, it turns out that the (g — 2)^ does not place any tight constraints on the 
models we consider here. Hence, all models we show in the later sections are compatible 
with (g - 2) M . 



3. Models for the dark matter halo 



The distribution of dark matter particles in the halo of the Galaxy plays a major role when 
making predictions for dark matter detection rates. Unfortunately, that distribution is 
rather poorly known, and one has to rely on large extrapolations. The approach we follow 
here is to take N-body simulations of hierarchical clustering in cold dark matter cosmologies 
as a guideline. Numerical results indicate that dark matter halos can be described fairly 
accurately by a universal density profile of given outer slope and singular towards the 
galactic center (there is a general consensus on the r -3 scaling at large radii, while there 
is still an on-going debate regarding how cuspy the profiles are), and by only two extra 
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parameters, e.g. a length scale and the normalization of the density profile at that scale, or 
better, in a language which is more appropriate in a structure formation study, the virial 
mass of the object M v i r and its concentration parameter c v i r (roughly speaking, the total 
mass and some measure of what fraction of this mass is confined in the inner portion of the 
halo). Actually, since the first systematic studies [31], it has emerged that M v i r and c 



■mr 



are, at a given redshift, highly correlated, with smaller objects being more concentrated, 
an effect which has been understood in terms of the form of the power spectrum of density 



perturbations and of the redshift of collapse of a given class of objects [32, 33 1. 

N-body simulation results provide snapshots of halos before the baryon infall; the 
formation of the luminous components of the Galaxy is likely to have induced a back- 
reaction on the structure of the dark halo as well. We will consider two opposite and at the 
same time extreme regimes for this effect. In the first, the baryon infall occurs as a smooth 
and slow process, with no net transfer of angular momentum between baryonic and non- 
baryonic terms. This is the regime in which the dark halo gets adiabatically contracted, a 
process in which, in the limit of spherical systems with unchanged local velocity distribution 
(as it happens, e.g., if all particles are placed on circular orbits), the mass distributions in 
the initial and final configurations are related by pM: 



Mi{n)n = \M b {r } ) + M C DM(r f )} r f . (3.1) 

Here, Mj(r), Mf,(r) and Mcdm(t) refer, respectively, to the mass profile of the halo be- 
fore the baryon infall (i.e. the form we assumed we can infer from N-body simulation 
results), the baryon component as observed in the Galaxy today, and the cold dark matter 
component in its nowadays configuration. 

In the opposite regime, a large transfer of angular momentum between the luminous 
and the dark components is assumed during the baryon infall. Mechanisms that induce 
a heating up of the population of cold dark matter particles sitting at the center of halos 
have been invoked to smooth out the inner CDM cusps, for which there is no evidence in 
real galaxies, or maybe even incompatibility, especially in the case of small dark matter- 
dominated objects, such as dwarf and low surface brightness galaxies []35| (but see also 
f36fl ). As there is no mechanism intrinsic in CDM to remove central cusps, in several 
analysis the focus has been moved to the possible role of baryons. In a recent and, to some 
extent, extreme model |J7]] baryons sink in the central part of halos after getting clumped 
into dense gas clouds: the halo density profile in the final configuration is found to be 
described by a profile with a large core radius of the type proposed by Burkert p£|]: 

P° 

(1 + r/a) (1 + {r/aY) 

a profile that has been shown to provide accurate fits of the rotation curves for a large 



sample of spiral galaxies [39]. 



The procedure we implement to select a halo model goes into several steps, which 



we summarize briefly here and are described in detail elsewhere [40 1 : To a given halo, 
labeled by a pair (M v i r ,c v i r ), with correlation between M v i r and c v i r of the form suggested 
by results in numerical simulations, we apply one of the two procedures for baryon infall 
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sketched above. The gas and stellar terms in the Galaxy are modeled as a superposition of 
components, i.e. a stellar bulge, bar and disk, plus two gaseous disks, whose morphologies 
are assumed in agreement with photometric studies in the Milky Way, and in analogy 
to the form observed in external galaxies. These components and the dark halo in its 
final configuration are required to match available dynamical information: among others, 
we apply constraints following from the motion of stars in the Sun's neighbourhood, total 
mass estimates following from the motion of the outer satellites, and we require consistency 
with the Milky Way rotation curve and with measures of the optical depth towards the 
galactic bulge. 

As a result we extract two sample dark matter halo models: i) a Burkert halo profile, 
with M v i r = 1.3 x 1O 12 M0 and c v i r = 16 (we use for the definition of c v i r the convention of 



Bullock et al. [32] with a slightly different normalization to the mean virial overdensity with 



respect to, e.g., Ref. [31]), and described by Eq. ( |3.2| ) with a = 11.7 kpc and a local halo 
density pb^o) = 0.34 GeV cm -3 ; ii) an adiabatically contracted profile, derived starting 



from a profile with the non-singular form extrapolated in Ref. [41] from simulations with 
the highest resolution so far, and M v i r = 1.8 x 10 12 Mq and c v i r = 12 (we label the final 
halo as N03 profile). The form of the latter is obtained numerically with Eq. ([O]) and we 
cannot give its explicit form: the local halo density in this case is /9/vo3( r o) = 0.38 GeV cm -3 
and the profile has a pronounced cusp towards the galactic center, driven in its most inner 
portion by the r -1 ' 8 singular scaling in the bulge profile 42], which dominates over the non- 
baryonic dark component. Adiabatic contraction makes the profile even steeper than the 
r _1 cusp in the profile originally proposed by Navarro, Frenk and White [31], and actually 



the result obtained starting from this or shallower profiles is equivalent. The approximation 
of adiabatic contraction is assumed to be valid up to the radius of 1 pc, which corresponds to 



the range within which the mass budget gets dominated by the central black hole [43]. We 



consider a scenario for the black hole formation in which the CDM system is perturbed and 



a core in the density profile of about 1 pc appears (see, e.g., [44], f45|), the opposite regime 



compared to the case of adiabatic growth of the black hole which induces the formation of 



a very dense spike within the innermost parsec 46]: the latter is essentially equivalent to 
adding a point dark matter source at the Galactic center, and can be treated separately. 

The final step is to derive a velocity distribution function consistent with the dark 
matter density profile considered. For isotropic velocity distributions and in the limit of 
spherically symmetric systems, the cases we restrain to in this analysis, the distribution 
function F is determined univocally by the density profile through Eddington's formula [47]: 



F(e) 



2 



d 2 pf (ffi 1 f dp 



(3.3) 



where VD^r) = — <3?(r) + <I>(r = oo), with $>(r) the potential for the Galaxy (i.e. including 
all components as determined from the dynamical information, in a toy model where the 
density profiles of stars and gas are assumed to be spherical as well, an approximation which 
may seem somehow drastic but has little influence on the final result), e = —E + §{r = 
oo) = —Ek in + *(r), and E and E^n, respectively, the total and kinetic energy. The 
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functions F(e) for our two sample models are computed numerically, after changing variable 



of integration in Eq. (3.3) from $ to the radius, r, of the spherical system. 

Details of the two sample models we consider here are give in Ref. p(J and tabulated 
forms of these density profile and distribution functions are also included in the latest 
release of the DarkSUSY package [12]. We just underline here that, for the first time, 
dark matter detection is treated in a scheme where dark matter halo profiles are directly 
compared to available dynamical constraints and fully self-consistent density profiles and 
local velocity distributions are considered, so that self-consistent comparisons between 
direct and indirect detection techniques are possible. 

4. Detection methods 

We will consider the discrimination and detection prospects with most of the techniques 
that have been proposed so far for WIMP dark matter searches. We describe briefly here 
the main underlying assumptions, and how we compare, in each class of experiments, with 
current sensitivities and with future ultimate-goal experiments. 

4.1 Direct detection 

Concerning direct detection, we will present our results in terms of the neutralino-nucleon 
spin independent scattering cross sections (more precisely, we will plot the neutralino- 
proton cross section, as the neutralino-neutron one is analogous). Spin dependent cross 
sections are not considered here, because the detection prospects are less encouraging. The 
reader should keep in mind that in the c x - p which is provided, is encoded, besides the de- 
pendence on SUSY parameters, a slight dependence on the underlying effective Lagrangian 
approach and on nucleonic matrix elements. Regarding the first issue we use approximate 
heavy squark expansions rather than full one-loop gluonic interactions. On the latter, we 
take a standard set of parameters ^| (e.g., for the s-quark operator, the most impor- 
tant in determining the scattering cross section, we set the parameter fj, = 0.14, slightly 



smaller than the value implemented in other analyses, see [12, |5(| for details). The scatter- 
ing cross sections are compared to the sensitivity level reached by current direct detection 
experiments and to the projected final reach of future large scale (one-ton-size) detectors 
(unfortunately, in the scheme we consider we do not find any model which produces cross 
sections at the level needed to provide an annual modulation effect compatible with the 
effect found by the DAMA Collaboration |5l|). For definiteness, we will refer to the EDEL- 
WEISS 2002 [|5^] exclusion limit, and to the sensitivity foreseen in the proposed XENON 
detector . Both of these are recalculated for the local velocity distribution functions in 
the two halo models considered here, taking into account relevant effects, such as target 
materials and form factors, thresholds, exposures and backgrounds (sensitivity curves for 
the standard Maxwell-Boltzmann velocity distribution have been derived as well, and cross 
checked against published results, with which they agree). 

4.2 Neutrino flux from the Sun and the Earth 

The neutrino flux induced by pair annihilations of neutralinos trapped at the center of the 
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Sun, is computed according to the standard procedure described in Ref. [54], except for the 



treatment of the neutralino capture rate, which is derived numerically implementing the 



expressions in [55| which are valid for a generic velocity distribution function (rather than 
the approximate expression appropriate for a Maxwell-Boltzmann velocity distribution Q , 
which are used by most authors), applied to the two cases considered here. This allows 
for a fully consistent comparison between this signal and direct detection. The results 
are presented in terms of muon fluxes, above the threshold of 1 GeV, generated by the 
muon neutrino flux, and compared for reference with the best exclusion curve among those 
provided by currently operating neutrino telescopes. As the current best limits are from 
the SUPER-KAMIOKANDE Collaboration in 2002 @ we will use these for the current 
sensitivity, and for the next generation of neutrino telescopes we will use the projected sen- 
sitivity of the km 2 -size detector which is being built by the ICECUBE Collaboration [|57j . 
For these projected sensitivities, we will use the expected sensitivity for the particular 
annihilation channel that dominates the neutralino annihilation, i.e. the sensitivities will 
be worse for soft channels (e.g. bb) than for hard ones (e.g. W + W~, Z°Z° and tt). The 
neutrino flux induced by neutralinos captured by the Earth has been computed as well. 
However, our results show that the expected signals are always orders of magnitude below 
the expected sensitivities. One should also remember that the Earth signal is highly cor- 
related to the size of the spin-independent coupling to nucleons (in the same way the Sun 
signal is correlated to the spin-dependent and spin- independent coupling), and we find that 
direct detection seems always more promising. We are not going to consider this signal 
further. 

4.3 Antimatter fluxes 

Neutralino pair annihilations in the Galactic halo may produce a significant amount of 
antimatter, even at the level of the standard secondary flux generated in the interaction 
of primary cosmic rays with the interstellar medium. We consider here neutralino induced 
fluxes of antiprotons, positrons and antideuterons. The prediction for these fluxes consists 
of several steps: our particle physics setup fixes the pair annihilation cross section at zero 
temperature <T ann t; (i.e. the probability, at the present time, for non-relativistic neutralinos 
to annihilate), and the branching ratios for the various annihilation channels. The frag- 
mentation and/or decay of these annihilation products give rise to the stable antimatter 



species. We have modeled this process with the Pythia [58] 6.154 Monte Carlo code, sim- 
ulating for each allowed two-body final state and for a set of 18 neutralino masses, except 
for D sources for which we have implemented the prescription suggested in Ref. [|59| to 
convert from the p-n yields. The strength of local sources scales with the number density 
of neutralino pairs locally in space, i.e., in terms of the dark matter density profile, with 
1/2 (p x (x ) / m x ) 2 , where m x is the neutralino mass. The choice of halo profile enters then 
critically: large enhancements in the source functions are provided by the adiabatically 
contracted profile which is cuspy towards the Galactic center. 

The next step is to model the propagation of charged cosmic rays through the Galactic 
magnetic fields. We consider a two-dimensional diffusion model, in which reacceleration 
effects are not included explicitly, but mimicked through a diffusion coefficient which takes 
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the form of a broken power law in rigidity, R, 



D = D (R/R ) os if R > R 

D = D if R<R . (4.1) 



E.g., in Ref. |60|| such a form has been tested with the Galprop |51| propagation code, 
showing that, with an appropriate choice of parameters, it is possible to reproduce fairly 
well the ratios of primary to secondary cosmic ray nuclei. We adopt here the same pa- 
rameter setup, namely Eq. ( |4,1| ) with Dq = 2.5 x 10 28 cm 2 s _1 and Rq = 4 GV, in a 
cylindrical diffusion region of radius equal to 30 kpc and half height equal to 4 kpc, plus a 
galactic wind term. In the case of antiprotons and antideuterons this is interfaced into the 
diffusive-convective code illustrated in Ref. |62| , which neglects energy losses (the particle 
is removed whenever a scattering with a nucleus takes place). For positrons, we instead use 
the code developed in Ref. [p3], improved and extended to allow for the implementation 
of a diffusion coefficient in the form of Eq. ( |4.1| ), and to keep a full two-dimensional struc- 
ture. In this code there is no convective term but there is a term accounting for positron 
energy losses due to inverse Compton scattering on starlight and the cosmic microwave 
background. 

The quality of the data on the local antiproton and positron cosmic ray flux has kept 
improving in recent years. We will compare the predicted fluxes to the antiproton data 
collected by the BESS experiment during its flights in 1997, 1998, 1999 and 2000 @ 
with fairly good statistics in the energy range between 180 MeV and 4.2 GeV, and by 
the CAPRICE experiment during its 1998 flight in the range between 3 and 50 GeV. 
For the positron fluxes, we refer to data published by the HEAT Collaboration on their 



1994-1995 flight §§, by the CAPRICE team regarding a flight in 1994 f67|, plus data 
obtained by MASS-91 in 1991: the overall energy range covered by these measurements 
extends from 460 MeV to 34.5 GeV. We have chosen not to include in our analysis data 
which have been reported just as antiproton or positron fractions (rather than absolute 
fluxes) and datasets such as the one on positrons from the AMS test flight |68| mapping a 
low energy interval in which a primary neutralino-induced contribution is not likely to be 
present. 

To compare with these data it is necessary to include the effect of solar modulation, 
i.e. the effect of propagating fluxes from interstellar space to our location inside the solar 
system and against the solar wind. To sketch this effect, we implement the one parameter 



model based on the analytical force- field approximation by Gleeson & Axford 6£] for a 
spherically symmetric model. The solar modulation parameter, sometimes dubbed Fisk 
parameter <S>p 7C], is for simplicity assumed to be charge-sign independent, and is derived 



fitting the measured local proton cosmic-ray flux to the solar activity corresponding to the 
period of data taking of each of the considered datasets. Finally, it is also necessary to 
consider a model for the secondary antimatter fluxes, which play the role of backgrounds in 
our analysis. For this task we refer again to the Galprop code, running the code in the 
same setup we use for the signals and for the current best estimate for proton cosmic-ray 
flux, and extracting the output on secondary antiprotons and positrons. For both species, 
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these backgrounds provide excellent fits of the data: we obtain, for background only, a 
reduced \ 2 equal to 0.82 for antiprotons and to 0.95 for positrons. 

Future measurements of antimatter in space, such as with the PAMELA satellite de- 
tector |7l]] and the AMS spectrometer [72] on board the International Space Station Alpha 



(ISSA), with instruments that will operate over very long exposures (tentatively, at least 
3 years), will provide data on an energy range which is wider than that of present data 
and with much better statistics. To address the perspectives of discrimination of the sig- 



nal against the background, we follow the approach outlined in Ref. [73] and define the 
quantity: 

h= L„ " e ^hw (4 - 2) 

where <& s (i£) and <&b(E) refer, respectively, to the signal and the background flux, and the 
integral extends over the energy interval in which the integrand is non negligible. This 
represents the continuum limit, up to an overall factor which accounts for the exposure 
times effective area of a given future experiment, of the \ 2 variable of the form: 

1 s+ b , oj (4.3) 

(AA*) 2 

under the assumption that, in each energy bin i, the number of signal events A] is sub dom- 
inant with respect to the number of background events A?, which in turn nearly matches 
the number of observed events A*, and that errors can be approximated as statistical er- 



rors, i.e. A A* ~ v A* ~ y Nfc. We have assumed as well that the measurement covers 
the energy range where the bulk of the signal is expected, and that the background is 
known, as one can (optimistically) expect from future determinations of the propagation 
parameters from, say, high precision measurements of ratios of secondaries to primaries 
for several light nuclei. With these hypotheses, we will indicate up to about what level 
of the parameter the PAMELA detector will be able to reject the presence of a signal 
contribution in the antiproton and positron flux. 

The case for antideuterons is different because in this case, restricting to a low energy 
window, the background flux is expected to be negligible [58], and even detection of 1 



event would imply discovery of an exotic component. Regarding the detection prospects, 
we will consider, as the ultimate reach for an experiment in the future, that of the gaseous 
antiparticle spectrometer (GAPS) ]74| . This is a proposal for an instrument looking for 
antideuterons in the energy interval 0.1-0.4 GeV per nucleon, with estimated sensitivity 
level of 2.6 x 10 _9 m _2 sr _1 GeV~ 1 s~ 1 , to be placed either on a satellite orbiting around the 
earth or on a probe to be sent into deep space. 

4.4 Gamma rays 

A certain amount of gamma-rays is also produced as a result of the fragmentation and 
decay of annihilation products from neutralino annihilations in the halo (at the same time a 
monochromatic component may be present, due to two-body annihilation states containing 
a photon; this component is not considered here, see, e.g., [75], for a discussion of detection 
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prospects). The induced flux is generally small, unless one deals with a dark matter halo 
profile with sharp density enhancements locally in space. As one of the two profiles we are 
considering here is cuspy towards the Galactic center (GC), we will compute the neutralino- 
induced 7-ray flux in that direction. 

The EGRET experiment, on the Compton Gamma-Ray Observatory has resolved a 7- 
ray source towards the GC |7^j , tentatively extended (~ 1.5°, of the order of the EGRET 
angular resolution) rather than point-like, and with a spectrum that is sensibly harder 
than the spectrum expected for the diffuse 7-ray flux due to the interaction of primary 
cosmic rays with the interstellar medium. Quite intriguingly, the measured flux shows the 



distortion that the diffuse flux would have in case of a WIMP-induced component [77]. 
We will compare the fluxes we find in our framework with the intensity and spectrum of 
the EGRET 7-ray source. Note, however, that alternative tentative explanations for this 
source have been proposed [fTSfl , and that even the identification of the location of the 
source with the GC has been questioned [79|. Should one of these two issues be confirmed 



by the upcoming observations with the GLAST satellite |8C], more stringent limits than 
those shown here could be derived. 

More recently, observations at higher energies of a 7-ray flux in the GC direction with 
ground-based Air Cerenkov Telescopes ^] have been reported. Given the neutralino 
mass range we will consider here, we find that the EGRET measurement gives always a 
more stringent bound. 

5. Results and discussion 

In line with most previous analyses, we will sample the 5-dimensional mSUGRA parameter 
space choosing a few values of tan/3 and Aq, and varying mi/ 2 an d mo for both sign(fi). 
We consider three different regimes: i) tuq < and Aq = 0; ii) small my 2 , moderate mo 
and large Aq; and iii) mo 3> m^ 2 . For each of these regimes we find the isolevel curves for 

X 



Q, x h 2 = 0.103, and select the parts of these curves that fulfill the accelerator constraints. 



5.1 Bulk, slepton coannihilation and funnel regions 

In the regime at mo < m^ 2 and Aq = there are three distinct frameworks which drive 
the relic density of binos to be in the cosmologically favored range. Starting at small m^j 
and small mo, we find the "bulk" region, where the relic density of pure binos is set by their 
annihilation strength into fermions mediated by a sfermion; it is a region, however, which 
is almost entirely excluded by accelerator constraints, see, e.g., the discussion in Ref. |jS3fl , 
Going to larger mi/2 there are two possibilities: for small or moderate values of tan/3 the 
curves with a fixed value of the relic abundance enter into the slepton coannihilation tail, 
a thin strip on the border with the region where a stau is the LSP and the relic abundance 
gets dominated by coannihilation effects with the lightest stau, as well as with the lightest 
selectron and smuon (this region was recognized in Ref. [Q| and the relic density calculation 



has been discussed in several recent analyses, including ]85|, |86|, |87], |88|, |89|, |Uj). As sample 
cases of this regime, we have searched for relic density isolevel curves for tan (3 = 10 and 
tan (3 = 30. As already mentioned, for large tan /3, one hits instead the region where the 
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Figure 1: The expected direct detection spin independent scattering cross section (left panel) and 
neutrino- induced muon fluxes from the Sun (right panel) for a few sample cases at mo < mi/2- 
Along each of the curves we have fixed fl x h 2 = 0.103. The lower curves for the direct detection 
limits are for the N03 profile, and the upper for the Burkert profile. The predicted muon fluxes are 
for the Burkert profile, with the N03 fluxes being only marginally larger. 



pair annihilation of neutralinos goes on the s-channel resonance of Hq and/or Hq: Q x h 2 
isolevel curves get two separate branches, at slightly larger and slightly smaller values of 
mo with respect to the value that, at a fixed m]y 2 , drives the neutralino mass to be one 
half of the Higgs mass; the two branches merge at the largest allowed values of mi/2 (and 
hence of the neutralino mass) for a given value of the relic density. This is more evident 
for sign(fi) = — 1, as in this case a sharp resonance appears in the neutralino LSP region 
even for values of tan (3 around 40 or so, while for positive /j, one has to go to very large 
tan/3. We consider the case tan/3 = 53, sign(fi) = +1 and tan/3 = 46, sign(fj.) = — 1; for 
both of them the isolevel curves for f2 x /i 2 in the cosmologically interesting band, start, at 
small rrii/2, in the bulk region, going to larger m^, enter the upper resonance branch, and 
then bend down into the resonance branch with lower mo, going now to smaller 7771/2, an d 
finally falling into the slepton coannihilation tail, as for the case with lower tan f3. Although 
we have derived, for all tan/3 and both sign{n), the isolevel curves in the full plane, we 
will have to take into account that the b — > S7 constraints will cut a large portion of the 
parameter space, especially at large tan/3 and negative fi. Also in the case of sign(fi) = — 1 
the SUSY contribution to (g — 2)^ is negative. It is still within our broad range of allowed 
values, but slightly disfavoured. 

Sitting on the isolevel curves for Q x h 2 = 0.103, we show our results in figures where 
we plot the neutralino mass on the horizontal axis, in the range allowed in the lower end 
by the accelerator constraints, the H® mass limit for positive /x and b — > S7 for negative 
//, and in the upper end by the fact that the isolevel curve hits the region where the stau 
is the LSP. The neutralino mass labels univocally the location on an isolevel curve, except 
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along a portion of the funnel stripe (as mentioned above). In Fig. [T] we show the expected 



spin independent scattering cross section and muon rates induced by the neutrino fluxes 
from the Sun. Regarding the first, we notice that as the squarks are heavy, the main 
contributions to the cross section always come from diagrams with t-channel exchange of 
H® and H® (which are almost decoupled from each other, and then, in their turn, coupled 
mostly, respectively, to u-type SM fermions and d-type SM fermions). The size of these 
two contributions are comparable. The one through H f rises with tan (3, as the H® d d 
vertex scales like yl + (tan (3) 2 . The two contributions have a constructive interference 
for positive /i, while they have a destructive interference for negative /x, giving very small 
values of the cross section. Even for positive \x and large tan f3, the scattering cross sections 
are rather small, at least a factor of 50 below the current sensitivities; this feature is due 
to the fact that we are dealing with rather pure bino-like neutralinos, while the neutralino 
coupling to the CP-even Higgs is roughly scaling with Z g (l — Z g ), with Z g being the gaugino 
fraction, i.e. Z g ~ \Nu\ in our case: the coupling is suppressed for models with Z g too 
close to 1. Considering the detection prospects for the future, for fi > there is a fair range 
of models which lie above the projected exclusion curve of the XENON experiment (the 
two dashed curves refer to the two halo model considered, with the lower one being the 
one for the N03 profile, mainly because the local halo density is larger in that case), but 
the high mass end, for each tan/3, will not be tested. For negative (jl, none of the allowed 
models will be testable. Regarding the neutrino fluxes from the Sun, we have compared 
with the current limit from Super-Kamiokande |5(| (derived for an average type of neutrino 
energy spectrum). We have also compared with the expected sensitivity for IceCube p7fl , 
derived for a soft annihilation channel as the dominating annihilation channel is bb for 
these models. As Fig. [I] shows, the rates are too low even for future large size neutrino 
telescopes. The flux displayed is for the Burkert profile, but even for the N03 case the 
increase is marginal. The fluxes are much lower than in more favorable scenarios, because 
the spin dependent cross section is lower than in these cases (and hence the capture rate 
is lower), as it is mediated, for pure binos, mainly by squark s-channel exchanges, and 
squarks in this scenario are rather heavy. 

The case for indirect signals from halo annihilation is introduced in Fig. ^. As a 
preliminary remark, we stress that we do not expect correlation of the halo rates with spin 
independent or dependent scattering cross sections, or better, we may find large rates for 
large cross sections (and actually, to some extent, we do) but it is never the same effect 
driving them to be large, i.e. crossing symmetry is never at work here. In the left hand side 
of Fig. 121 we plot the annihilation rate at zero temperature divided by the neutralino mass 
squared, which is the quantity appearing in the source function of all halo rates. Naively, 
implementing an approximate solution of the Boltzmann equation describing neutralino 
decoupling in the early Universe in which temperature dependences in the annihilation 
cross section are neglected, a ann v scales with the inverse of the relic abundance H: 



We are working at a fixed r2 x /i 2 = 0.103; the corresponding curve is shown as a green 
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Figure 2: In the left panel we show the annihilation cross section divided by the neutralino mass 
squared for a few sample cases in the limit mo < miw. In the right panel we compare the predicted 
antiproton and gamma-ray fluxes against current data, see the text for details. 



dotted line in the figure. If this approximation holds, the curves for the models we are 
considering should align on top of it, while we see that this is marginally true just for the 
isolevel curve following a trajectory within resonances (which are very broad and hence 
thermal effects are not dramatic). For the other curves two effects come in: first there is 
the S-wave suppression of the cross section for the annihilation of non-relativistic Major ana 
fermions into light fermions, an effect that appears already for small neutralino masses; 
on top of this, going to larger masses, the thermally averaged annihilation cross section 
setting the relic abundance gets dominated by coannihilation effects. In this case, we can 
see that the scaling with mass of (cr ann v) / 'm^ on the curves referring to the tan/3 = 10 or 
30 cases is much more rapid than 1/rrvt. In the figures are plotted also curves which give 
a feeling for the constraints from current data. They are calculated under the hypothesis 
that the bb annihilation channel dominates, a very good approximation, especially for the 
funnel case. The dashed lines correspond to approximate 3a exclusion curves based on 
current antiproton data, for the two halo model considered, the lower curve corresponding 
to the N03 profile (and assuming the propagation model we picked is the correct one); 
they are nearly flat because we have factored out most of the dependence on mass in the 
signal in the quantity we plot on the vertical axis. The dotted lines are the analogous for 
positrons, and we see that, in this case, the curves are well above the predictions even for 
the N03 profile. Finally, the dash-dotted lines indicate cross sections for which the EGRET 
excess from the GC is better fitted by a neutralino-induced component in case of the N03 
profile plus a (normalization free) diffuse background component (see [77] for details on the 
background and on how this fit is done); the case for the Burkert profile does not fit into 
the range of cross sections displayed on the vertical scale. In the lower mass tail of the case 
with tan j3 = 53 and the N03 profile, it looks that current data may put a constraint on the 
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Figure 3: Perspectives to detect dark matter within the mSUGRA framework and the too < m i/2 
regime. Dotted curves in the upper part of the two figures refer to approximate 3 a exclusion curves; 
in the lower parts of the figures, models with visibility ratios larger than 1 should be considered 
detectable in future experiments. Models with il x h 2 — 0.103 and two different values of tan/3 are 
displayed. 



model. This is shown on the upper right-hand side of Fig. ^, where we plot the reduced x 2 
for the antiproton flux with respect to current data, and the ratio between the predicted 
gamma-ray flux in the energy bin between 4 GeV and 10 GeV and the flux measured by 
EGRET in the same bin (which is the highest one for the EGRET source, and turns out to 
be the one setting the tightest constraint). For the N03 profile the lowest mass range seems 
to be excluded by both signals; note however, that the 3a and la exclusion lines we plot 
for the antiprotons do not take into account uncertainties in the propagation models, and 
should then probably be slightly relaxed (but most likely not up to a reduced x 2 OI 3), and 
that the estimate of the gamma-ray flux is very sensitive (much more than the antiproton 
flux) to details in the halo profile towards the GC: for the Burkert profile, the flux dies out 
to 4 orders of magnitude below the limit (for reference, in the language of Ref. frEfl , (J), 
the line of sight integral angularly averaged over a cone with a 1.5° aperture, is equal to 
about 7 x 10 4 for the N03 profile considered here). 

In Fig. ||, we address the detection prospects for future experiments, for the case 
tan/3 = 53 and tan/3 = 30, both with positive /i. In the lower part of the figures are 
plotted the ratio between the spin independent scattering cross section and the projected 
sensitivity of the XENON proposed experiment, for the two halo models, as well as the 
predicted D flux, at the earth position and in a period towards solar maximum, over the 
sensitivity of GAPS (on a satellite close to Earth). As already mentioned, direct detection 
looks rather promising, and also very promising seems, for the N03 profile, but partially 
also for the Burkert profile, the measurement of the D flux (which would be even more 
promising should the experiment be performed with a detector placed on a probe sent 
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Figure 4: The direct detection scattering cross section (left panel) and the expected neutrino- 
induced muon fluxes from the Sun (right panel) for our chosen stop coannihilation region. The 
curves are analogous to those in Fig. [I]. 

into deep space). In the upper part of the figures, with the same mass scale, we plot the 
quantity 7$ defined in Eq. (|4.3|), in units of cm~ 2 s _1 sr -1 , compared with the 3 a exclusion 



curves in this quantity for an experiment with an effective area of the size of PAMELA [90] 
(see for details) and an exposure of 1 yr or 3 yr. As can be seen, these indirect channels 
looks less promising, and the positron channel seems always worse than the antiproton one 
(but it should be kept in mind that we are assuming that the background is known; if a 
signal is indeed present, it might be easier to spot it in the positron flux at relatively high 
energy, rather than in the antiproton flux at low energy). 

5.2 Stop coannihilation region 

This is the regime at small moderate mo and large |Aq|. We focus here on a sample 

case in which the parameter A§ is tuned to generate a thin slice where neutralino-stop 
coannihilations drive the relic abundance into the cosmologically interesting range: we 
choose Aq = —3000 GeV, tan (3 = 5 and sign(/i) < 0. The stop coannihilation strip opens 
up at small m\/2 and significantly larger mo, in a region where the neutralino is still a 
rather pure bino and has a mass of about 100-300 GeV. In many cases, the LEP constraint 
on the H® boson mass severly affects these regions and can rule out the entire region with 
relic densities in our desired range. For the parameters we have chosen here, the Higgs 
constraint comes in at lower masses and sets a lower limit on the neutralino mass of 208 
GeV. 

In Fig. | we show the direct detection scattering cross section and the expected muon 
fluxes for these models together with current and expected future sensitivities. The dom- 
inating annihilation channel for these models is gluon-gluon below the ti threshold and ti 
above, i.e. a soft spectrum below m t and a hard above. We have used this to plot the 
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Figure 5: In the left panel we show the annihilation cross section divided by the neutralino mass 
squared (red solid line) for our chosen stop coannihilation region. Also shown are the current limits 
coming from p, e + and 7 measurements. The lower set of these curves are for the N03 halo profile 
and the upper for the Burkert profile. In the right panel, we show the reduced \ 2 f° r the fits to 
antiproton data for these models and the visibility of gamma rays from annihilation at the galactic 
center. 




Figure 6: Future detection prospects for our chosen stop coannihilation region. The lower set of 
curves are for the N03 profile and the upper for the Burkert profile. The direct detection scattering 
cross section is too low to be seen even by future detectors. 



expected sensitivity for IceCube |??J (the jump at mt just reflects the change from soft to 
hard annihilation spectrum). As is evident from these figures, none of these models can be 
detected with neutrino telescopes and direct detection experiments. 
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In Fig. | we instead focus on signals from annihilations in the galactic halo. In this 
regime, neutralino pair annihilations are, above threshold, dominated by the ti channel, 
which is not S-wave suppressed and it is a rather copious source of antimatter particles. 
When ti is not kinematically allowed, below the ti threshold, the gluon-gluon channel 
dominates, but this is never the case for the models displayed as the set of models below 
the threshold and beyond are excluded by the mass constraint from LEP. Annihilation 
to gluon-gluon below the threshold is though assumed in the exclusion curves shown. As 
can be seen in the left panel of the figure, we are now much closer to the current sensitivities, 
but still not quite to a level of excluding models. In the right panel, we check explicitly 
the corresponding reduced x 2 f° r antiproton fluxes against current data, which falls within 
the l-o" range for almost all models. We also compare with the gamma ray flux from the 
galactic center, showing that the predicted fluxes are always lower than the flux measured 
by EGRET for the GC source, even for the cuspy N03 profile. 

We now move on to future detectors, and in Fig. || we show the detection prospect 
by future searches. The antiprotons are marginally detectable with future detectors, like 
PAMELA in our cuspy N03 profile, whereas the positron fluxes are too low. The an- 
tideuteron fluxes, on the other hand, are high enough to be detectable by e.g. the GAPS 
probe, especially if the halo profile is cuspy, like the N03. 

To summarize, the stop coannihilation region seems to be best probed by the an- 
tideuteron fluxes as measured by e.g. a future GAPS probe. The direct detection and 
neutrino telescope rates are far too low to be seen even by anticipated future detectors. 
The reason that annihilation in the halo is more advantageous, is probably that annihilation 
into ti dominates and this gives rather large fluxes of cosmic rays, whereas the scattering 
rates are low since the crossed diagram is not significant for the scattering cross section. 

5.3 Focus point region 

We describe now the regime at ttiq S> rn,\/2- Sitting on a isolevel curve at fixed relic 
abundance in this regime essentially means that, at each mass we are selecting a preferred 
Higgsino fraction for the lightest neutralino. At a fixed Higgsino fraction, the annihilation 
cross section into gauge bosons, which dominates for pure Higgsinos, decreases as the mass 
increases, hence the higher in mass we go, the more pure Higgsinos we need to have to 
compensate for this decrease. In this respect, tan/3 do not play much of a role, and the 
cases for positive and negative \i are perfectly analogous (except that again the (jl < 
case is slightly disfavored by (g — 2)). We choose to show, as sample cases, slices of the 
parameter space at = 0, fi > and tan/3 = 30 or tan/3 = 50. 

The predicted direct detection spin independent scattering cross section and neutrino- 
induced muon fluxes from the Sun are shown in Fig. ^. The spin independent scattering 
cross section is in this case dominated by one single diagram, the one mediated by H®, as 
the squarks are extremely heavy, and H® turns out to be very heavy as well: it follows 
that, as can be seen, there is not much of a dependence on the tan/3 parameter, contrary 
to the cases at small mo we have discussed. Also, the cross sections are (relatively) large 
over most of the mass range, as there the Higgsino-Bino mixing in the lightest neutralino 
is fairly large, and we already pointed out a large mixing enhances the coupling of the 
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Figure 7: The expected direct detection spin independent scattering cross section (left panel) 
and neutrino-induced muon fluxes from the Sun (right panel) for two sample cases in the limit 
mo » mi/2- Along each of the curves we have fixed Q x h 2 — 0.103. 



neutralino to CP-even Higgses. Just at a mass of around 1 TeV we get to the limit of 
very pure Higgsinos and the cross section at that point drops dramatically. Except for 
this corner, this framework looks testable with the XENON detector (while unfortunately 
even this case is unconstrained by current experiments). For the expected sensitivity in 
IceCube, we have used the fact that for these models the neutralinos annihilate mainly 
to Z°Z°, W + W~ or tt, i.e. they have hard annihilation spectra. The signal due to the 
neutrino flux from the Sun looks rather promising as well in a future perspective: the effect 
which is driving it to levels which may be tested at future large size telescopes is the fact 
that the spin dependent scattering cross section is rather large, driven by the diagram with 
a Z° boson exchange (with Higgsino Z° coupling that does not depend on tan /?, but again 
just on Higgsino fraction and mixing between the two Higgsino states, hence with the two 
cases at different tan/? nearly overlapping). The two halo profiles change slightly the level 
of the flux, with about a 23% increase for the N03 profile, an effect which is mostly due to 
the higher local density (which increases accordingly the capture rate), but with the two 
different velocity distribution functions playing some minor role. 

As shown in Fig. ||, the annihilation cross section of neutralinos in halos today is rather 
accurately set by the relic abundance, with curves at different tan (3 overlapping again and 
nearly aligning over the curve from which a ann v is estimated using Eq. flB.ip . Comparison 
with current data sets is performed in exactly the same way as for the other cases we 
discussed. Except for the antiproton flux in case of low masses and for the N03 profiles, 
current data do not set constraints on the configurations we have chosen. 

The prospects for detection of these dark matter candidates in the future are shown in 
Fig. ^. The focus point region, and in particular the slice of it we have selected at fixed value 
of the relic abundance and therefore with large Bino-Higgsino mixing, is the region which 
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Figure 8: In the left panel we show the annihilation cross section divided by the neutralino mass 
squared for two sample cases in the focus point region. In the right panel we compare the predicted 
antiproton and gamma-ray fluxes against current data, see the text for details. In the latter we 
display just the case for the N03 profile, as for the Burkert profile signals are much smaller than 
measured fluxes. 




Figure 9: Future detection prospect in the focus point region, for the tan/3 = 50 case we discussed 
in the text. The tan/3 = 30 case is perfectly analogous. 

looks more promising from the point of view of dark matter detection. Direct detection will 
test models up to very large masses, cross checks with very clean signatures may come as 
well from neutrino telescopes, and the measurement of D flux (for small neutralino masses, 
even in case of the Burkert profile) . If the halo is cuspy, the neutralino- induced antiproton 
and positron flux may be at a detectable level as well, and might be soon singled out in 
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upcoming experiments. 



6. Conclusion 

We have performed a detailed analysis of current limits and detection prospects of neu- 
tralino dark matter in the mSUGRA framework. We have focused on models with a thermal 
relic density, as estimated with the DarkSUSY numerical package, in the currently favored 
cosmological range, and considered all relevant regimes in the parameter space. Direct and 
indirect detection rates have been computed implementing two dark matter halos, with 
fully consistent density profiles and velocity distribution functions, and opposite histories 
for the transition between the stage of a CDM halo prior to the baryon infall and a halo 
embedded in a galaxy with inner portion dominated by the luminous components, as is 
the case for the Milky Way halo. This has allowed, for the first time, a fully consistent 
comparison between direct and indirect detection. 

In general, we can conclude that most of the mSUGRA models considered here are 
not excluded by any of the current dark matter searches. For some models (low mass stau 
coannihilation region and low mass focus point region), we overproduce antiprotons and 
gamma rays from the galactic center in our cuspy N03 profile (but not with the cored 
Burkert profile). 

For future experiments, we have found that in the region of small mo, direct detection is 
rather promising if [i is positive and tan (3 is large, a feature due to the scattering amplitudes 
mediated by CP-even Higgs bosons summing coherently and to the coupling in the d d 
vertex becoming large. In the same region, but for different reasons, the neutralino-induced 
antiproton, positron and especially antideuteron fluxes could be detectable. In the stop 
coannihilation region, both the direct detection and the neutrino telescope rates are too 
low to be detectable even with future experiments. The most promising technique to test 
these models is to search for an antideuteron flux with an experiment like GAPS; large 
fluxes follow in this case from large annihilation rates into top quarks. Finally, in the funnel 
region, direct detection looks very promising because of the large portion of both Bino and 
Higgsino in the lightest neutralino. An eventual signal in direct detection experiments may 
be cross checked with the measurement of the induced neutrino flux from the Sun, and 
may even be anticipated through measurements of cosmic ray antimatter fluxes; both of 
these kinds of signals are expected to be large because of the large Higgsino fraction in this 
region. 

As we have stressed, this analysis applies to one specific framework, and relies on 
specific effects emerging in that framework. Extrapolations to other perfectly viable models 
are possible, once the corresponding relevant effects are singled out in those models as well. 
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